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Rhenium diboride is a recently recognized ultra-incompressible superhard material. Here we 
report the electronic (e), phonon (p), e-p coupling and thermal properties of ReE>2 from hrst- 
principles density-functional theory (DFT) calculations and neutron scattering measurements. Our 
calculated elastic constants (cii = 641 GPa, C12 = 159 GPa, en = 128 GPa, C33 = 1037 GPa, and 
C44 = 271 GPa), bulk modulus (B w 350 GPa) and hardness (H w 46 GPa) are in good agreement 
with the reported experimental data. The calculated phonon density of states (DOS) agrees very 
well with our neutron vibrational spectroscopy result. Electronic and phonon analysis indicates 
that the strong covalent B-B and Re-B bonding is the main reason for the super incompressibility 
and hardness of ReB2. The thermal expansion coefficients, calculated within the quasi-harmonic 
approximation and measured by neutron powder diffraction, are found to be nearly isotropic in a 
and c directions and only slightly larger than that of diamond in terms of magnitude. The excellent 
agreement found between calculations and experimental measurements indicate that first-principles 
calculations capture the main interactions in this class of superhard materials, and thus can be used 
to search, predict, and design new materials with desired properties. 

PACS numbers: 71.20.-b, 62.20.Dc, 63.20.-e, 65.40.-b 
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Hard materials are of great scientific interest due 
to their numerous technological applications. Unfortu- 
nately, almost all superhard materials (diamond, cubic 
BN etc.) are expensive because they either occur nat- 
urally in limited supplies or have to be made at high 
pressure synthetically. Therefore, intense research efforts 
have been carried out to design superhard materials 1 . 
Recently, it was found that rhenium diboride can be syn- 
thesized at ambient pressure with potentially low cost, 
and the resulting ReB2 crystal has super incompressibil- 
ity along the c axis, comparable to that of diamond, and 
high hardness, comparable to that of cubic BN. 2 The 
mechanical properties of ReB2 were also correctly pre- 
dicted by a recent theoretical work 3 -. To more fully un- 
derstand this unusual material, we have used a combined 
first-principles and neutron scattering study to further 
investigate the electronic, elastic, phonon and thermal 
properties of ReB 2 . 

Our calculations were performed within the plane- wave 
implementation of the generalized gradient approxima- 
tion (GGA) to density-functional theory (DFT) in the 
PWscf package^. We used Vanderbilt-type ultrasoft po- 
tential with Perdew-Burke-Ernzcrhof exchange correla- 
tion. A cutoff energy of 680 eV and a 16x16x6 fc-point 
mesh (generated using the Monkhosrt-Pack scheme) were 
found to be enough for total energy to converge within 
0.01 meV/atom. Spin-polarized calculations resulted in 
zero spontaneous spin polarization for the material in- 
vestigated; thus here we focus on the results from spin- 
restricted calculations. 

Neutron scattering measurements were undertaken at 
the NIST Center for Neutron Research. The inelastic 
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FIG. 1: (Color online) (a) A unit cell of ReB 2 . (b) The 
electron localization function of the (1 1 -2 0) crystal plane. 
The four types of chemical bonds are denoted using dashed 
lines. The strong directional, covalent B-B (labeled "1") and 
Re-B (labeled "2" and "3") bonding are apparent. 



neutron scattering (INS) measurements were performed 
with the BT-4 filter-analyzer neutron spectrometer. The 
neutron powder diffraction was performed with the BT- 
1 high-resolution powder diffractometer. ReB2 powder 
sample used in both experiments was synthesized using 
the method of direct heating of the elements (Re and 
n B) in vacuum at 1273 K, as reported in Ref.2. The 
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purpose of using n B isotope is to avoid the large neutron 
absorption of normal boronS. 

RcB2 has a simple hexagonal structure (space group 
PQs/mmc) with experimental a — 2.900 A and c = 7.478 
AS, as shown in Fig. 1(a). The two Re atoms occupy the 
sites (1/3, 2/3, 1/4) and (2/3, 1/3, 3/4) while the four 
boron atoms occupy the sites (2/3, 1/3, z), (1/3, 2/3, 
1/2+z), (2/3, 1/3, 1/2-z), (1/3, 2/3, l-z) with experi- 
mental z — 0.048. We first optimized both the lattice 
parameters and the atomic positions of the ReB2 struc- 
ture. The relaxed parameters (without considering the 
zero point motion), a = 2.9007 A, c = 7.4777 A, and z 
= 0.0478, all agree very well with the experimental val- 
ues. We note that the crystal structure of ReB 2 is quite 
different from other metal diboride compounds, a ma- 
jority of which assume MgB 2 type structure. In MgB 2 , 
boron atoms are on the same plane, with a graphcnc 
structure^, while they are buckled in ReB2. In addition, 
Mg is located above the center of B-hexagon in MgB 2 , 
while in ReB 2 , Re is located right on top of the B atom. 
The structural configuration of ReB 2 is mainly due to 
the strong hybridization of Ke-d and B-p orbitals, which 
we will discuss in detail below. 



TABLE I: The calculated bond population, bond length, and 
real-space force constants for selected pairs of ions in a (1,1,- 
2,0) plane in ReB2. The positive (negative) bond popula- 
tion indicates bonding (antibonding) character. Note that the 
force constants are very large for nearest neighbor B-B and 
Re-B pairs, indicating very strong local covalent bonding. 
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No. 


Bond 


Population 


Length (A) 


Force Constant (eV/A x ) 


1 


B-B 


1.96 


1.820 


4.55 


2 


Re-B 


0.72 


2.257 


4.34 


3 


Re-B 


-0.25 


2.227 


5.08 


4 


B-B 


-0.14 


3.025 


1.36 



To understand the high hardness of ReB 2 , it is es- 
sential to look at its electronic structure. According 
to the calculated valence charge distribution, there ex- 
ist large electron densities between two neighboring B 
atoms, and between the Re atoms and its neighboring B 
atoms, indicating strong directional, covalent B-B and 
Re-B bonding. This is as expected, since a material 
with large hardness must contain highly directional, short 
and strong bondsS*^. To get more detailed informa- 
tion about the bonding nature, we performed electron 
localization function (ELF) analysis 10 and Mulliken bond 
population analysis 11 . In particular, the ELF plot is very 
useful in terms of distinguishing different bonding inter- 
actions. The value of ELF is in the range of to 1 by 
definition. High ELF value corresponds to a low Pauli ki- 
netic energy, as can be found in covalent bond. A value of 
ELF near 0.5 corresponds to delocalized electron density 
as found in metallic bonding. The ELF plot of the (11-2 
0) crystal plane of ReB 2 is shown in Fig. 1(b). The dom- 
inant feature is the rather strong B-B bond (labeled "1" ). 
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FIG. 2: (Color online) (a) The calculated electronic band 
structure of ReB2, along high-symmetry directions in the Bril- 
louin zone. T, A, H, K, M and A represent k points (0 0), (0 
0.5), (-1/3, 2/3, 0.5), (-1/3, 2/3, 0), (0, 0.5, 0) and (0, 0.5, 
0.5) respectively, (b) The total and partial density of states 
(DOS) of ReB2. The main feature of the electronic structure 
of ReB2 is a strong hybridization of Re-d and B-p orbitals, 
which results in the strong covalent Re-B bond. 



The two Re-B bonds (labeled "2" and "3") clearly have 
different bonding characters. The B-B bond along the c 
axis (labeled "4") has a large length, but ELF shows that 
its contribution to the overall bonding is not negligible. 
More quantitative bond population analysis (see Table I) 
confirms that the B-B bonding is very strong and highly 
covalent, with essentially a double bond nature. Inter- 
estingly, the Re and its two neighbored B atoms sitting 
right above and below it (along the c axis) form anti- 
bonding, while Re forms strong covalent bonding with 
all other six B neighbors. In table I, the real-space force 
constants between these ions, obtained from the phonon 
calculation (discussed later), are also shown, which indi- 
cate strong B-B and Re-B bonding as well. The B-B and 
Re-B bonds, all together, form a covalently bonded three 
dimensional network. 

In Fig. 2(a), we show the electronic band structure of 
ReB 2 . The total density of states (DOS) and the partial 
DOS projected onto atomic orbitals are shown in Fig. 
2(b). Clearly, the material is metallic and the electronic 
structure of ReB 2 near the Fermi level is governed by a 
strong hybridization of Re-e? and B-p orbitals, which re- 
sults in the strong covalent bonding of B-B and Re-B. 
In terms of this bonding characteristic, ReB 2 is similar 
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FIG. 3: Least-square fit of the total energy data vs. strain 
parameters. The dots are from the first-principles calculations 
and the solid lines are the least-square fit. The deformation 
matrices for each distortion are also shown. The five hexag- 
onal elastic constants of ReB2, along with the bulk modulus 
and shear modulus are summarized in the right-bottom panel. 



to the previously well studied OsB2r^^J^J^ although 
the latter has an orthorhombic structure. We also note 
that the density of states at the Fermi level, N(E^) = 
1.6 states/eV, is quite high and slightly larger than that 
of the 40 K superconducting MgB2 8 , raising the question 
of possible superconductivity in ReB2, which we will ad- 
dress later. 

We next discuss the elastic constants of the relaxed 
ReB 2 structure, which was calculated directly. For 
hexagonal crystal, there are five independent elastic coef- 
ficients: en, C12, C13, C33 and C44. Therefore, we applied 
five symmetry-independent strains 1 -: (7, 7, 0, 0, 0, 0), 
(7, -7, 0, 0, 0, 0), (0, 0, 7, 0, 0, 0), (0, 0, 0, 0, 2 7 , 0), (7, 
7, 7, 0, 0, 0) to extract the five unknowns. We used 7 
of ± 0.005, ± 0.0075, ± 0.01, ± 0.015. From the least- 
square fit of the total energy vs. strain datail (Fig. 3), 
we found that c n = 641 GPa, c 12 = 159 GPa, c 13 = 128 
GPa, c 33 = 1037 GPa and c 44 = 271 GPa. The quality 
of the data fit is excellent. By doubling the strain mag- 
nitudes, the fitted elastic constants only change slightly, 
suggesting that the anharmonicity is negligible in this 
structure. Note that our calculated elastic constant along 
c axis, C33 = 1037 GPa, is very large and indeed compa- 



rable to that of diamond (w 1040 GPa, the largest known 
elastic coefficient 18 ). This result conforms well with the 
reported experimental observation, namely, that the in- 
compressibility along the c axis is equal in magnitude 
to the linear incompressibility of diamond 2 -. Also note 
the anisotropy in the compressibility of the two differ- 
ent lattice directions, i.e., a axis is substantially more 
compressible than c axis, which is also consistent with 
experimental observation. For a hexagonal crystal, the 
Voigt bulk modulus B = 2(cn + c i2 + 2ci 3 + c 33 /2) / 
9, thus we obtain B — 350 GPa, in good agreement with 
the experimental value 2 of m 360 GPa. The Voigt shear 
modulus G = (7cn - 5ci2 - 4ci3 + 2C33 + I2C44 ) / 9, 
from which we obtain G — 283 GPa. Comparing our 
calculated elastic coefficients with Hao et. al.'s compu- 
tational results^, we found only a small difference (< 5 
%), which is likely due to the slight difference between 
different plane wave DFT codes. 

Materials with high bulk and shear modulus are of- 
ten hard materials. The hardness of a covalent or ionic 
compound can be directly calculated using the method 
recently proposed by Simnek and Vackai^, based on an- 
alyzing bond strengths and their densities. As discussed 
earlier, in the ReB2 crystals, each Re atom forms eight 
bonds with B atoms, six of which have dji e —B — 

2.257 A 

(bonding) while other two have slightly shorter d^e-B = 
2.227 A (anti-bonding). In addition, each B atom forms 
three covalent bonds with its neighboring B atoms (ds-B 
= 1.820 A). The Re- Re interaction has a weak metallic 
nature with a nearest neighbor distance of 2.90 A, thus 
their contribution to the hardness is negligible. The ref- 
erence energies^ (number of valence electron divided by 
the radius which makes the atom neutral) for Re and B 
were calculated to be 4.878 and 3.09 respectively. Using 
equation (6) in Ref.19, we derived H = 46.0 GPa. Al- 
though the anisotropy of the hardness was not considered 
here, the calculated average value falls right within the 
range of the experimentally measured hardness (« 30 to 
56 GPa)2. 

Next, we discuss the phonon structure of ReB2, which 
also reflects the mechanical properties. In addition, the 
lattice dynamics determine a wide range of other macro- 
scopic behaviors such as thermal and transport proper- 
ties, and the interaction with radiation (e.g., infra-red 
absorption, Raman scattering or inelastic neutron scat- 
tering). We performed the dynamical calculations on 
the optimized ReB 2 structure using the supercell method 
with finite differenced A 2x2x1 supercell was used, and 
the full dynamical matrix was obtained from a total of 
8 symmetry-independent atomic displacements (0.02 A). 
The unit cell of ReB 2 contains six atoms, which give rise 
to a total of 18 phonon branches. The phonon modes at 
r are classified as 



r(g = 0) = A lg + 2A 2u (IR) + 2B H 



B 



2u. 



E 2u 



2E 2g (R) + 2E lu (IR) + E lg (R), 



where R and IR correspond to Raman- and infrared- 
active, respectively. The crystal symmetry implies six 
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Raman- and six 7i?-active modes. In Table II, we list the 
calculated energies at T and corresponding mode charac- 
ters. 

The computed phonon dispersion curves and phonon 
density of states are shown in Fig. 4(a). The acoustic 
branches have steep slopes, indicating large elastic coeffi- 
cients. Using the low energy part of the acoustic branches 
(i.e., "elastic limit"), we can estimate the sound velocity 
(V) and thus also elastic constants. For example, for 
a hexagonal crystal, c 33 = p V (longitudinal [0001]) 2 and 
c 44 = pV (transverse [0001]) 2 , where p is the mass den- 
sity. We derived c 33 = 948 GPa and c 44 = 278 GPa. The 
slope determination close to q = certainly has some er- 
ror bar; nevertheless, the above estimated values are still 
in reasonable agreement with the more precise numbers 
calculated earlier by directly applying strains. 

The optical phonon branches are clearly divided into 
two groups. The three lowest energy optical modes (< 30 
meV) are dominated by Re motion while the high energy 
optical modes (> 50 meV) are dominated by the lighter B 
atom displacement. Note that for similar types of atom 
motion, we found that those along c axis always have 
higher energies than the one within the basal a — b plane. 
This is due to the covalent Re-B bonding along the c 
axis. For Boron motion within the plane, those involving 
B-B bond stretching (i.e., the motion of the bonded B-B 
pair being out-of-phase) are of higher energy than the in- 
phase motion, due to the strong B-B covalent bonding. 
Hence, the phonon structure is totally consistent with 
our ealier electronic structure analysis. 

TABLE II: The calculated phonon energies E and the char- 
acteristics of the optical modes of ReB2 at F-point. Note that 
in the table, "in-phase" and "out-of-phase" refer to the move- 
ment of the bonded B-B pair in the unit cell. (IR) and (R) 
indicates the IR and Raman active modes, respectively. The 
last column indicates the electron-phonon coupling for each 
modes at F-point, which are quite small compared to those of 
MgB 2 . 



Species 


E (meV) 


Dominant character/type 




2E 2g 


18.6 


Re in a-b plane 





Big 


28.5 


Re along c 


0.007 


2E 2u 


50.1 


B in a-b plane, in-phase 





2-Ei-u 


59.6 


B in a-b plane, in-phase (IR) 


0.008 




78.1 


B along c, in-phase (IR) 





2E\g 


85.2 


B in a-b plane, out-of-phase (R) 


0.002 


B 2u 


87.6 


B along c, out-of-phase 


0.025 


2E 2g 


90.4 


B in a-b plane, out-of-phase (R) 


0.016 


Big 


91.1 


B along c, out-of-phase 


0.021 




97.7 


B along c, out-of-phase 


0.018 



The calculated phonon structure is further validated by 
the neutron spectroscopy measurement. The measured 
INS spectrum is essentially the phonon DOS weighted 
by the neutron scattering cross sections of the elements. 
As shown in Fig. 4(b), the agreement between the exper- 
imental data and the calculated spectrum based on the 
incoherent approximation 2 ^ is excellent. 
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FIG. 4: (a) Calculated phonon dispersion curves along high- 
symmetry directions in the Brillouin zone for ReB2 and the 
phonon density of states, (b) Inelastic neutron scattering 
spectrum from ReB2 (at 4 K) along with the calculated spec- 
trum. The agreement is excellent. 

It is interesting to note that the phonon structure 
of ReB2 is quite similar to that of MgB2, a high-T c 
material^. With a MgB2-like phonon spectrum and a 
slightly higher N(Ei?) than that of MgB2 as mentioned 
earlier, one may wonder if RcB2 can exhibit superconduc- 
tivity at a "high" temperature. Experimentally, in the 
Re-B system, Re 3 B and ReyB 3 were found to have a T c of 
4.8 K and 3.3 K, respectively* 2 ^ We are not aware of any 
experimental measurement of T c in ReB2. We thus cal- 
culated the electron-phonon coupling parameters at V for 
ReB2 and tried to estimate T c . The results are shown in 
Table II. The coupling mainly takes place for B phonons 
along c axis. Compared to MgB2, the electron-phonon 
coupling in ReB2 is much weaker. Using the McMil- 
lian expression 2 ^ and only considering the zone center 
phonon, T c of ReB2 is estimated to be less than 1 K. It 
would be interesting to confirm this estimate experimen- 
tally. 

Finally, we discuss the thermal expansion of the ReB2 
structure, a property important for practical applica- 
tions. The temperature dependence of the lattice param- 
eters a and c are calculated within the quasi-harmonic 
approximation. The results were obtained by minimiz- 
ing the Helmholtz free energy 

F(a,c,T) = V(a,c) + £^{^(q) 

j q 

+ kTln(l-e- h ^^/ kT )\, 
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• c, exp. 
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FIG. 5: (Color online) (a) The calculated (lines) and ex- 
perimental (dots) temperature dependence of the structural 
parameters a and c of ReB2. (b) The calculated (lines) and 
experimental (dots) linear thermal expansion coefficient of 
ReB2. For comparison purpose, the calculated linear thermal 
expansion coefficient of diamond is also shown. The thermal 
expansion of ReB2 is nearly isotropic and the magnitude is in 
the same order as that of diamond. 



where the first term is the ground state energy and the 
second term is obtained by summing the phonon modes 
over the wavevectors in the Brillouin zone. In the quasi- 
harmonic approximation, the effect of the anharmonicity 
in the lattice energy is treated by allowing the phonon 
frequencies to depend on the lattice parameters. Hence, 
for a given temperature T, we first take a and c, mini- 
mize the atomic positions, and then calculate the phonon 
spectrum, and the free energy. Repeating this for other 
values of a and c, we find the optimum values of the lat- 
tice parameters that minimize the free energy at a given 
temperature. In this way, one obtains the temperature 



dependence of a and c. 

The calculated temperature dependence is shown in 
Fig. 5(a), along with the experimental lattice constants 
measured by neutron powder diffraction. Although the 
DFT-GGA calculation predicts slightly larger lattice con- 
stants than the measured values, the overall agreement 
on the lattice expansion is reasonably good. We found 
that the linear thermal expansion coefficients along a and 
c axes are both very small [see Fig. 5(b)], w 6.5 xlO -6 
K" 1 at 1000 K, slightly larger than that of diamond (« 
4.5xl0~ 6 KT 1 at 1000 K). Additionally, a and c expan- 
sion coefficient curves nearly overlap with each other, 
suggesting that the material is highly isotropic in terms 
of thermal expansion, which is preferred in most appli- 
cations. The excellent combination of hardness, incom- 
pressibility and small thermal expansion coefficient in- 
dicates that ReB2 has great potential for many appli- 
cations, such as material for cutting tools and abrasion 
resistance coating. 

In summary, we have applied density-functional theory 
and neutron scattering techniques to elucidate the elec- 
tronic, elastic, phonon, and thermal properties of ReB2. 
Our calculated elastic constants, bulk modulus and hard- 
ness are in very good agreement with the experimental 
data. Our electronic and phonon results confirmed that 
the strong covalent B-B bonding and Re-B bonding play 
a critical role in the incompressibility and hardness of 
ReB2. Our calculations indicate that ReB2 has a very 
similar phonon spectrum to MgB2 with a comparable 
N(Ei?). However, we found a very small electron-phonon 
coupling, which suggests a very modest superconduct- 
ing temperature. The thermal expansion coefficient is 
found only slightly larger than that of diamond. The 
combined excellent mechanical and thermal properties 
suggest great potential for ReB2 to be used as a cutting or 
coating material. The excellent agreement found between 
DFT calculations and experimental measurements indi- 
cate that first-principles calculations are able to capture 
the main interactions in this class of superhard materi- 
als, and thus can be used to search, predict, and design 
other materials with properties (hardness etc.) better 
than diamond. 
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